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Abstract 

The quantum mechanics status of the probabihty vector current density has long seemed to be marginal. 
On one hand no systematic prescription for its construction is provided, and the special examples of 
it that are obtained for particular types of Hamiltonian operator could conceivably be attributed to 
happenstance. On the other hand this concept's key physical interpretation as local average particle 
flux, which flows from the equation of continuity that it is supposed to satisfy in conjunction with the 
probability scalar density, has been claimed to breach the uncertainty principle. Given the dispiriting 
impact of that claim, we straightaway point out that the subtle directional nature of the uncertainty 
principle makes it consistent with the measurement of local average particle flux. We next focus on the fact 
that the unique closed-form linear-superposition quantization of any classical Hamiltonian function yields 
in tandem the corresponding unique linear-superposition closed-form divergence of the probability vector 
current density. Because the probability vector current density is linked to the quantum physics only 
through the occurrence of its divergence in the equation of continuity, it is theoretically most appropriate 
to construct this vector field exclusively from its divergence — analysis of the best-known "textbook" 
special example of a probability vector current density shows that it is thus constructed. That special 
example in fact leads to the physically interesting "Ehrenfest subclass" of probability vector current 
densities, which are closely related to their classical peers. 

Introduction 

The quantum mechanical probability vector current density concept has long been at best hazily understood. 
Although special examples of probability vector current density have been obtained for particular types of 
Hamiltonian operator [1], there is no systematic prescription for constructing it, so its special examples could 
conceivably be attributed to happenstance. Its essential feature is supposed to be that, in consequence of the 
quantum mechanical conservation of probability, it jointly with the probability scalar density satisfies the 
equation of continuity, which compellingly suggests that it physically represents local average particle flux. 
It has, however, been claimed that the measurement of local average particle flux breaches the uncertainty 
principle because such a measurement involves arbitrary particle localization while it simultaneously yields 
information concerning particle velocity, and for that reason probability vector current density cannot have 
its only truly natural physical interpretation, namely that of local average particle flux [1]. This contention, 
if valid, would raise two linked vexing issues: (1) if probability vector current density cannot in principle be 
physically interpreted as local average particle flux, then what is its correct physical interpretation, and (2) 
is that correct physical interpretation still something that is of actual interest to physicists? 

Very forturnately, however, while it is indeed the case that measurement of local average particle flux 
restricts particle position while delivering information regarding particle velocity, these things do not in fact 
occur in such a manner as to challenge the uncertainty principle. To measure one of the vector components 
of local average particle flux at a given point, one passes a plane which is perpendicular to that vector 
component through that point. One then selects an arbitrarily small region of that plane centered on 
that point and measures the average rate that particles pass through that selected planar region. This 
average rate, divided by the area of the selected planar region, is an approximation to that particular vector 
component of average particle flux at that particular point, an approximation which is, in principle, refined 
by additionally shrinking the planar region on which that point is centered. The flux one thus measures 
does indeed reflect the component of average particle velocity in the direction perpendicular to the plane, 
provided, of course, that one knows the average particle density. The uncertainty principle, however, does 
not restrict the accuracy with which a particular component of particle momentum can be determined due 
to restrictions that are imposed on the particle's position in directions which are perpendicular to that 
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momentum component. To see this one need look no further than the Dirac canonical commutation rule, 



[xi,Xj]=0, \pi,Pj]=0, [xi,pj]=ih5ij, i,j = 



1,2,3 



and recall that the uncertainty principle only applies to quantized dynamical variable pairs that fail to com- 
mute, which is clearly not the case for a position component and a momentum component that are mutually 
perpendicular. Therefore the measurement of local average particle flux , notwithstanding that it restricts par- 
ticle position while delivering averaged particle velocity information — namely in the direction perpendicular 
to the plane in which particle position is restricted does not challenge the uncertainty principle. 

Thus there is in fact no known valid reason to refrain from physically interpreting probability vector 
current density as local average particle flux, even as the equation of continuity ostensibly satisfled by 
probability vector current density in conjunction with probability scalar density so cogently suggests is sound. 
The bona fide issue here is a different one, namely the lack of a general prescription for constructing that 
desired probability vector current density which actually satisfies the equation of continuity in conjunction 
with the probability scalar density. 

We therefore in the rest of this article concentrate on systematically working out such a probability vector 
current density for any canonical Hamiltonian operator whatsoever (at least in principle — as one can well 
imagine, anything worked out to such a degree of generality yields formidable expressions that cannot be 
expected to have transparency as their strong suit). We begin with the divergence of the probability vector 
current density, which, from the equation of continuity, must be equal to the negative of the time derivative 
of the well-defined probability scalar density. That time derivative, in turn, is the state-vector expectation 
value of {i/Ti) times the commutator of the set of position projection operators with the Hamiltonian oper- 
ator. Because there exists a linear superposition technique for the unique closed- form quantization of any 
classical Hamiltonian function [2] , it turns out that for Hamiltonian operators which have classical Hamil- 
tonian function antecedents we can also uniquely reduce the just-mentioned expression for the divergence of 
the probability vector current density, which involves the state-vector expectation value of a Hamiltonian- 
operator commutator, to a linear-superposition closed form. Thus the divergence of the probability vector 
current density can always be fully and uniquely worked out when the Hamiltonian operator has a classical 
Hamiltonian function antecedent, i.e., when it is a canonical Hamiltonian operator. 

The probability vector current density itself of course is naturally mathematically amhiguous, e.g., it 
tolerates the addition of an arbitrary vector field that is a curl. However, since its divergence is a state- 
vector expectation value (of a set of operators), it is physically entirely reasonable to restrict the probability 
vector current density to being a homogeneously linear functional of that divergence so that it likewise will be 
a state-vector expectation value. Indeed, the larger idea of restricting the probability vector current density 
to the very barest minimum that is compatible with its divergence makes impeccable sense because it is only 
that divergence which makes actual contact with the system's quantum physics, doing so, of course, through 
the equation of continuity that it satisfies in conjunction with the probability scalar density. Therefore the 
probability vector current density is devoid of quantum physical information that is not already implicit in 
its divergence. We shall formally implement this crucial point by stipulating not only that the probability 
current density is homogeneously linear in its divergence, but that it furthermore must not depend on any 
constants which are additional to those that are intrinsic to its divergence, and also that among the forms 
which are mathematical candidates for this vector field, only the most symmetric are to be considered because 
these add no further information to that available from its scalar divergence. These three stipulations result 
in a unique closed-form expression for the probability vector current density in terms of its unique divergence, 
and we shall see that the "textbook" best-known special example of a probability vector current density [1] 
is indeed consistent with these three stipulations. 

We now turn to the presentation of the ingredients that are needed to work out the probability vector 
current density in the manner that has just been outlined. As is well-known, the conservation of probability 
in quantum mechanics [1] follows from the Hermitian property of the Hamiltonian operator, namely ^ {t) = 
H{t), in conjunction with the Schrodinger equation. 



ihd\tl^{t))/dt = H{t)\xp{t)), 



(la) 



which together imply that, 



-iM{i){t)\/dt = {i){t)\H{t) 



(lb) 
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and therefore that probabiUty is conserved, 



dm)m))idt = m) \ [{i/n)H{t) + {-i/n)H{t)\ m)) = o. (ic) 

In the case of canonical quantum mechanics, H{t) is uniquely obtained from a classical Hamiltonian 
function Hci{[x, p],t) by the imposition of the following self-consistent slight extension of Dirac's canonical 
commutation rule [2], 



[/i (x) + 51 (p) , /2 (x) + 52 (p)] = i?i '{/i(x)+5i(p)r/2(x)+52(p)}" , 



(2) 



where the { , } are the classical Poisson brackets, and the overbrace is used to denote quantization — the 

overbrace is only used here for that purpose because the traditional hat accent for denoting quantization is 
not sufficiently extensible. The Eq. (2) slight extension of Dirac's canonical commutation rule turns out to 
self-consistently completely resolve the Dirac operator-ordering ambiguity in favor of Born-Jordan operator 

ordering [2] , which is exactly the same operator ordering that is implicit in the Hamiltonian path integral [3] . 
The unique Hamiltonian operator which follows from the classical Hamiltonian function Hci{[x, p], t) via the 

application of Eq. (2) is, of course, denoted as ifci([x, p], , or, when explicit reference to the underlying 
classical Hamiltonian function p], isn't necessary, as simply the canonical Hamiltonian operator 

H{[%p],t). 

When such a canonical H{[x.,p],t) describes the quantum mechanics, probability conservation, given 
by Eq. (Ic), can be expressed in terms of the probability scalar density {{^p{t)\x) {x\ip{t))) because of the 
following expansion in the complete set of position states |x), 

d{^{t)\^{t))/dt = (d/dt) J (V'(t)|x)(x|i/'(i))d^x = 0, (3a) 

which relation might plausibly be expected to imply that the time derivative d{{tp{t)\x) {x.\tp{t))) / dt of the 
probability scalar density {{^p{t)\x.) {x\tp{t))) is a perfect differential in the vector variable x, namely that 
there exists a probability vector current density j(x, [\ip{t)),H{[x,p],t)]) which satisfies, 

d{{m\^){^\m))/dt + Vx • j(x, t; Mt)),H{[Sc, p],t)]) = 0. (3b) 

Eq. (3b) is the equation of continuity that the probability vector current density j(x, t; [\ip(t)), ([x, p], t)]) 
is supposed to satisfy in conjunction with the probability scalar density ((V'(f)|x)(x|i/;(i))) as a plausible 
consequence of the conservation of probability that is given by Eq. (3a). In light of the Schrodinger equations 
given by Eqs. (la) and (lb), as they apply to the particular case that H{t) = H{[x,p\,t), we can see that 
actually fulfilling the Eq. (3b) equation of continuity requires that, 



Vx-j(x,t;[|VXt)),i/([x,p],t)]) = -rf((^(t)|x)(x|^(t)))M = 

{i/h) [m)\x){x\Hi[%p],t)\m) - m)\H{[%p],t)\x){xm)} 
?fi\{2i/hmt)\^){^\Hi[^,p],t)m)) 



(3c) 



where the symbol 5i that occurs after the last equal sign of Eq. (3c) denotes the real part of the bracketed 

expression which follows it. We thus see that to obtain the divergence of j{x,t;[\tp{t)),H{[Si,p],t)]) wc 
need to evaluate the core part {x\H {[x.,p],t)\tp{t)) of that bracketed expression for an arbitrary canonical 
Hamiltonian operator i?([x, p],t). To achieve this goal it will be extremely useful to emulate the procedure 
of Ref. [2] wherein such an arbitrary H{[x,p],t) is linearly decomposed into "Fourier component" operators 
of the form g^'^^ '^e^^P '/'* because, of course, 

^x|e=F'^-ke±'P-i/fi|^(t)) = e=F'^-'*(x|e±'P-'/^|V'(t)) = e=F*''-''(x ± I|V'(t))- (3d) 



3 



The probability vector current density divergence by lineeir superposition 

Following Ref. [2], we note that the key step for decomposing H{[x, p],t) into "Fourier component" operators 
is the orthodox corresponding Fourier decomposition of its underlying classical Hamiltontonian function 
ifci([x, p], f). This, of course, follows from the identity, 

ffei([x, p], = / d^x' / d^p' 5(^|;x - x') 5(^^p - p') H,i([x', p'], i), (4a) 

after we insert into it the Fourier delta-function representations, 

5^^\-K - x') = (27r)-^ /d^k e-'k-(x-x')^ siN^-^ _ p') = (211%)-^ J d'^l e^Hp-p')/n_ (4b) 

Since 7Jci([x, p],t) is a real- valued function, after inserting Eq. (4b) into Eq. (4a) it is convenient to further- 
more explicitly discard the vanishing imaginary part of the result, which yields, 



ifci([x,p],t) = (47r2;i)-^/(i^x' /dVifci([x',p'],t) Jd^kjd^l X 
(cos(-x' • k p' • l/n) cos(-x • k -I- p • l/h) + sin(-x' • k p' • l/h) sin(-x • k -|- p • l/h)) . 

The fact that quantization is a linear process [2] permits us to conclude from Eq. (4c) that. 



(4c) 



H{[%p],t) = Fci([x,p],0 = (47r2;i)-^/d^x' /dVJ?ci([x',p'],0 Jd^k Jrf^l X 

/ ' ^ ' ^\ 

( cos(-x' • k -I- p' • l/h) cos(-x • k -I- p • l/h) + sin(-x' • k p' • l/h) sin(-x • k -|- p • l/h) j . 



(4d) 



As in Ref. [2], we obtain exp(— ix • k -|- ip • l/?i) from Eq. (2). From this we further immediately obtain 

, " s , " s , " ^ 

exp(ix • k — ip • l/h) , and those two results together yield cos(— x • k + p • l/h) and sin(— x • k + p • l/h) , 
which is what we require for insertion into Eq. (4d). We shall, however, insure that these cosine and sine 
quantizations are expressed as linear combinations of products of exponential operators of the "Fourier 
component" type that were introduced in Eq. (3d), because it is those "Fourier component" type operators 
which are transparently useful for the evaluation of the core part of the probability vector current density's 
divergence Vx • j(x, [\ilj{t)),H{\x,p],t)]), as is apparent from Eqs. (3c) and (3d). 

We now carry out the Eq. (2) quantization of exp(— zx • k + zp • l/7i) with the goal of expressing its result 
in terms of Eq. (3d) "Fourier component" type operators. With the needed exponential ingredients inserted, 
Eq. (2) reads, 

which, written out, is, 

g-ix-kgip-l/fi _ gip-l/fig-ix-k _ ^j-j, . j-j g-ix-k+ip-l/R _ (4f) 

Eq. (4f) can then be reexpressed as, 



^gi(k-l)/2 _ g-i(k-l)/2^ g-ix-k+i5-l/R ^ ^^jj. . g-ix-k+ip-l/fi _ (4g) 

which yields the unique exponential quantization, 



g-ix.k+ip.i/n ^ g-ix.k+i^.i/n sin(ik ■ l)/(ik • 1) = e-*^-''e'?-V'^e-*5k-i sin(ik • l)/(ik • 1), ('*^) 

where the form farthest to the right in Eq. (4h) is expressed in terms of the desired "Fourier component" 
type operator. We now make the simple substitutions k — >• — k and 1 — >■ —1 in Eq. (4h), which produce the 
additional useful result that. 



gix.k-ip-i/ft =e'^-k-i9-i/fisin(ik-l)/(ik-l) = e'''-''e-'P-V^e-'5k-isin(ik-l)/(i^ (^i) 
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Combining Eqs. (4h) and (4i) then yields, 

'cos(-x-k + p-l/;i)' = (1/2) (^e^i^-kgip-i/ft + eiJ-kg-i9-i/fi"j e-^i^-^sm{^k ■ l)/(ik • 1), (4j) 

and, 

sin(-x • k + p • l/h) ^ (l/(2i)) (^e-'^-'^e'P-'/^ - e^-^e-^P'^/^^ e-^i^'^ sin{^k ■ l)/(ik • 1), (4k) 

which are expressed in terms of the desired "Fourier component" type operators. We can now insert Eqs. (4j) 
and (4k) into Eq. (4d), which in turn is inserted into Eq. (3c), after which Eq. (3d) is straightforwardly 
applied. The upshot is the general linear superposition expression for the divergence of the probability 
vector current density that results from an arbitrary canonical Hamiltonian operator iJ([x, p],f) whose 
classical antecedent Hamiltonian function is ilci([x, p],t), 

Vx-j(x,t;[|V(t)),i?([S,p],t)]) = 
(47r2;i)-^/d^x' /d^p'(i7ci([x',p'],t)/?i) jd^k Jd^l (sin(ik-l)/(ik-l)) x 
(cos(-x' ■ k + p' ■ l/h) 5R[ie-*5k-i(^(t)|x) {e-''''-^{x + l\xP{t)) + e^^-'^(x - l\rp{t)))] + (^a) 

sin(-x' • k + p' ■ l/h) ^_e-'i^■\'^|J{t)\x) (e-'^-''(x + l\tlj{t)) - e*^-'^(x - 1|V'(0))] ) • 

One interesting special case of Eq. (5a) occurs when _ffci([x, p], f) has no dependence on the configuration 
variables x. In that case the integrations over the variables x' and k that occur in Eq. (4d) are obviously 
superfluous, and, if actually carried out in Eq. (5a), simply eliminate x' along with a factor of (27r)~''^, while 
setting k to 0. Thus in the case that i?ci([x, p],t) has no dependence on x, Eq. (5a) simplifies to, 

Vx • j(x, t; mt)),Hm,t)]) = (27r?i)-^ / d^p' iH,ii[p'],t)/h) J d^l x 

(cos(p' • l/h) mm\^) ((x + i\m) + (x - i\m))] + (sb) 
sin(p' . l/h) 3?[(v(t)ix) ((x + i\^{t)) - (x - imm) • 

Similarly, if iJci([x, p], t) has no dependence on the momentum variables p, then the integrations over p' 
and 1 in Eq. (5a) eliminate p' along with a factor of {2Trh)~^ , while setting 1 to 0. Thus in the case that 
Hc\{[x,p],t) has no dependence on p, Eq. (5a) yields, 

Vx-j(x,t; mt)),H{[9:],t)]) = {2n)-'' J d^x' {H,,{[x'],t)/h) J d^k x 

(cos(-x' • k) SR[jKV(t)|x)|2 (e-'"-'' + e'"-'')] + (5c) 
sin(-x'-k)3?[|(V'(0|x)|^ (e-i^-k_e»''-k)]^ =0, 

namely that Vx •j(x, t; [\^p{t)),H{[x.],t)]) vanishes identically. This can also be verified directly from Eq. (3c) 
for such a canonical Hamiltonian operator H{[x],t), because (x| is one of its eigenvectors, with the corre- 
sponding real eigenvalue Hci{[x.],t), and {ip{t)\x){x\t{j{t)) is real-valued as well. We therefore see that any 
terms of Hci{[x,p],t) which have no dependence on p can simply be truncated from Hci{[x,p],t) before 
Hci{[x, p],t) is inserted into the general linear superposition expression of Eq. (5a) for the divergence of the 
probability vector current density Vx ■ j(x,t; [\il;{t)) , H{[x,p],t)]). 

Our ultimate goal, of course, is the probability vector current density j(x, i; [\tp{t)), H{[x,p],t)]) itself, 
rather than its divergence. We now turn our focus to obtaining it, recalling that in the Introduction we set 
out three stipulations to be imposed on it to make it formally consistent with the fact that it is devoid of 
physical information that is not already implicit in its divergence. 

Obtaining the probability vector current density from its divergence 

A homogeneously linear form in terms of its divergence for the nth component of the probability vector 
current density is given by, 

{jix,t;mt)),H{[x,p],t)]))^ = 
Wn ( Vx ■i{xi,...,x'^,...,XN,t; [\il}(t)),H{^,p],t)\)dx'^, n = l,...,N, ^^^^ 
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where the weights Wn satisiy wi H Vwn = 1- We of course stipulated not only that the probability vector 

current density is homogeneously linear in its divergence, but also that it has no dependence on constants 
which are additional to those that are intrinsic to its divergence. The dependence of the expression on the 
right-hand side of Eq. (6a) on the set of constants {x^^^ , • • • , } that are its lower limits of integration is 
readily removed by setting all those lower limits of integration to — oo. That the result of thus introducing 
an infinite integration interval is well-defined (i.e., does not diverge) is directly tied in with the conservation 
of probability, as one can see from Eqs. (3a) and (3b) (with the latter repeated in Eq. (3c)). Finally, we 
have also stipulated that among the forms of the mathematical candidates for the probability vector current 
density, only the most symmetric are to he considered. That stipulation, and to a certain extent also the 
injunction against additional constants not intrinsic to its divergence, requires us to set every weight value 
Wn of Eq. (6a) to l/N, n= 1,. . .,N. With these specifications, Eq. (6a) becomes, 

(j(x,i; mt)),Hi[^,p],t)]))n = 

- (6b) 
{1/N) V.-j(a;i,..., a;;v,i;[|V'(*)),^^"([x,p],i)])rf<, n = l,...,N. ^ ' 

J —oo 

Eqs. (6b) and (5a) together are the linear superposition construction of the probability vector current density 
j(x, t; [|?/'(t)), i?([x, p], t)]) for any canonical Hamiltonian operator H{[x.,p],t) whose classical Hamiltonian 
function antecedent is i?ci([x, p], t) — see Eq. (5a) for the way that i?ci([x, p], t) is utiUzcd. 

We now turn to the best-known Hamiltonian operator, namely |p — po(x, i)P/(2m), for which direct 
application of Eq. (3c) suffices to extract the corresponding probability vector current density divergence 
Vx • j(x, i; [|i/)(t)), |p — po(x, t)p/(2TO)]) [1], which we abbreviate as Vx • j for convenience. It turns out 
that this divergence can simply be algebraically manipulated into the explicit form of a divergence operator 
Vx- acting on a certain vector field [1], without actually making use of Eq. (6b). We can then inquire 
whether the resulting "textbook" probability vector current density j(x, t; [\ip{t)), |p — po(x, i)p/(2TO)]) [1] 
that is specifically thus obtained by algebraic manipulation of its divergence Vx • j in fact adheres to the 
three stipulations used to derive Eq. (6b), namely homogeneous linearity in Vx • j, no additional constants 
beyond the ones intrinsic to Vx • j, and the maximum possible symmetry. 



Do the three postulated stipulations hold for the best-known special example? 

The key to applying Eq. (3c) to the Hamiltonian operator |p — po(x, i)p/(2m) is obviously to work out its 
Eq. (3c) core part, (x||p - po(x, i)p/(2m)|r/'(t)), 

(x||p - po(X,i)|V(2m)I^W> = (x|(|pP - P • Po(x,f) - Po(X,t) • P + \po{%t)\^)m))/{2m) = 

(x|(|p|2+iW;^.po(x,t) -2po(x,i) •p + |po(x,t)|2)|V'(i))/(2m) = (7a) 
[-^2 V2^(x, t) + 2i?ipo(x, t) ■ VxV(x, t) + {ihV^ ■ po(x, t) + |po(x, i)|2)^(x, t)] /(2m), 

where ipixJ-) '= (x|^i'(t)), and of course there is also its complex conjugate ^/'(x, t) =^ (^/;(t)|x). We now 
obtain the divergence Vx • j(x, t: [lipit]), |p — po(x, t)\'^ / {2m,)]), which we abbreviate as Vx • j for convenience, 
by merely putting Eq. (7a) into Eq. (3c). But we also find that with considerable additional algebraic 
manipulation that result can be explicitly presented as a divergence operator Vx- acting on a certain vector 
field, 

Vx • j = [i2i/n) |x) (x| |p - po (x, t) I V(2m) |^(i))] = 
m-iSR[-inv;(x,t)(V2V^(x,t)) -2V;(x,t)(po(x,0 • Vx^/'(x,i))- |^(x,t)|2(Vx-Po(x,i))] = 

(2m)-i[V;(x,i)((-»?i)V2V;(x,i)) + {{th)yii;{^,t))i^{^,t)]- 
m-i[V^(x,i)(po(x,t) • VxV(x,t)) + (po(x,t) • VxV^(x,t))^(x,i) + |^(x,t)|2(Vx • Po(x,t)] = (7b) 
(2m)-iVx • i)((-i?i)VxV'(x, t)) + iiih)V^^Pix, t))7/.(x, t)]- 
m-i[(Po(x,i) • Vx|V'(x,t)n + |^(x,t)|2(Vx • Po(x,t))] = 
Vx • [Vi(x,t)((-in)VxV(x,i)) + ((i7i)VxV^(x,t))^(x,t) - 2po(x,i)|^(x,t)|2]/(2m). 

Here the divergence operator Vx- has been simply factored out of the divergence expression Vx ■ j which 
follows from Eq. (3c). Therefore it is abundantly clear that the resulting "textbook" probability vector 



6 



current density [1], 



j(x,t;[|V(t)),|p-po(X,t)lV(2rn)]) = f7c) 
[VS(x,t)((-i;i)VxV(x,t)) + ((a)Vx^(x,i))V'(x,i) - 2po(x,i)|V(x,t)|2]/(2m), 

is homogeneously linear in its Eq. (3c) divergence Vx • j and also that it has no additional constants beyond 
those that are intrinsic to Vx • j- Inspection of its Eq. (7c) expression also reveals this particular probability 
vector current density to be highly symmetric; indeed its N-dimensional form is the completely symmetric 
generalization of its one- dimensional case. Thus this "textbook" best-known special example of a probability 
vector current density [1] definitely adheres to the three stipulations which underlie Eq. (6b). 

One naturally wonders which subclass of the class of canonical Hamiltonian operators H{[x.,p],t) has 
members which all emulate |p — po(x, t)|^/(2m) insofar as having the property that the divergence operator 
Vx- can simply be explicitly algebraically factored out of Vx ■ i{x,t:[\'(jj{t)), H{[x,p\,t)]). An important 
clue for working out that subclass is the fact that this property is manifest in the classical limit, where 
the single-particle vector current density can readily be shown through precisely such explicit algebraic 
factorization to be equal to the singular classical single-particle scalar density times the classical particle 
velocity. Interestingly, the straightforward quantization of the singular classical single-particle scalar density 
turns out to be equal to the quantum particle-position projection operator in the Heisenberg picture, and, 
of course, the quantum probability scalar density is sim,ply an expectation value of that particle-position 
projection operator. Thus it isn't greatly surprising that there exists a subclass of the canonical Hamiltonian 
operators H{[x,p],t) for which, in Ehrenfest-theorem style, the probability vector current density turns out 
to be equal to the expectation value of the quantization of the just-mentioned result for the classical single- 
particle vector current density, namely the singular classical single-particle scalar density times the classical 
particle velocity, a result that, exactly as in the classical case, is arrived at for this "Ehrenfest subclass" of 
the canonical Hamiltonian operators by the explicit algebraic factorization of the divergence operator Vx- 
out of Vx • j(x, t; Mt)),H{[Si, p],t)]). 

The classical vector current density and an Ehrenfest-like theorem 

As is familiar from its electromagnetic application [4], the classical single-particle scalar density is given by 
the singular expression, 

p,i{y,t) = 5^^U,i{t)-y), (8a) 
which, irrespective of the value of t, satisfies, 

Jpci{y,t)d^y=l, (8b) 

in precise analogy to the probability-conservation property of the quantum single-particle probability scalar 
density {tp{t)\y){y\tp{t)), which satisfies, 

j{m\y){ym))d''y = (mm)} = i- (sc) 

Therefore it is plausible that there exists a classical single-particle vector current density jci(yjt) which 
satisfies the equation of continuity in conjunction with pci{y,t), 

dPciiy, t)/dt + Vy • jci(y, t) = 0. (8d) 

For Eq. (8d) to hold it must be the case that, 

Vy ■M{y,t) = -dp,i{y,t)/dt = -{S(^\x,i{t)-y), i/ci([xci(t), Pci{t)],t)} = 

- ( Vx., it) S^^'l^cl {t) - y)) • Vp^, (t) //el ( [Xcl {t) , Pel (t)] ,t) = 

{Vy5(^\K,,{t) - y)) • Vp^,(,)i?ei([xci(i), Pcm,t) = ^^^^ 
Vy • [<5(^l;xci(t) - y) Vp^,(t)/fci([xci(t), Pci(i)],i)] = Vy ■ [S^^'lMt) - y) dx,i{t)/dt] , 

where { , } denotes the classical Poisson bracket, and we have used Hamilton's first classical equation of 
motion, dxci{t)/dt = Vp^j(()i/ci([xci(t), Pc\{t)],t). Eq. (8e) clearly shows that the divergence operator Vy- 



7 



explicitly algebraically factors out of the divergence of the classical single-particle vector current density 
Vy • jci(y, i), thus yielding for the classical single-particle vector current density jci(y,t) itself, 

^PciW-^ci([xci(t), Pc\{t)],t) — Pci{y,t) dxci{t)/dt. (8f) 

Prom Eq. (8f) we see that the classical single-particle vector current density is equal to the classical single- 
particle scalar density times the classical particle velocity, which is a simple, physically graphic result that 
turns out to have significant resonance in canonical quantum mechanics as well. 

The quantum/classical linkage arises both from the Eq. (2) relation between commutator and Poisson 
brackets and from the fact that the quantization of the classical single-particle scalar density Pdiyit) = 
5^^\^\{t) — y) is the Heisenberg-picture version of the quantum position projection operator |y)(y|, whose 
expectation value in the state is the quantum probability scalar density (V'(f)|y)(y|V'(t))- This aspect 

of the quantum position projection operator |y)(y| can be obtained from the fact that its application to any 
position eigenstate vector |x) yields, 

|y)(y|x) = 5(^\y - x)|y) = 5^^\^ - y)|x), (9a) 

which is identical to what the application of 5^^\Si — y) to that position eigenstate vector |x) yields, namely, 

5(^^x-y)|x) =<5(^^;x-y)|x), (9b) 

and the fact that, in the canonical quantum regime, the set {|x)} of position eigenstate vectors is complete. 
Therefore, the quantum probability scalar density {'tp{t)\y){y\tlj{t)) has the alternate expression, 

{m\y){y\m) = imis^'^l^ - y)\m), (qc) 

which is easily shown to produce the following alternate form of Eq. (3c), 

vyj{y,t;[\m),H{[%p],t)]) = {i/nmmsn^-y), m%p],mm)- (9d) 

Eq. (9d), in turn, implies the following alternate form of the combination of Eqs. (6b) and (5a), 

(j(y,i; [|VW),i?([x,p],i)]))„ = 

(1/^) (V?i)(V(OI[^^^lx - (yi, . . . , 2/;, . . . , yr,)), H{[x, p],t)]\m)dy'n, n = l,...,N. 

J — OO 



(9e) 



Eq. (9e) does not, however, provide information on the consequences of undertaking Pourier decomposition 
of the canonical Hamiltonian operator _ff([x, p],i), as Eq. (5a) docs. 

Much less general than Eq. (9c), but doubtless of greater physical interest, is the probability vector current 
density result for the subclass of canonical Hamiltonian operators H{[%p\,t) which satisfy a particular 
"Ehrenfest" requirement, namely that. 



[5^^\Si - y), ^([x, p], t)] = in {(5(^{x - y), H,,{[k, p],t)} , (l^^) 

where, of course, 

^([x,p],i) = Xi([x,p],i)'. W 
Prom Eq. (2) it is immediately apparent that canonical Hamiltonian operators of the form, 



H{[%p],t)= X([p],t) + y([x],t) =K{[p],t) + V{[Si],t), (10c) 

do indeed satisfy the "Ehrenfest" requirement of Eq. (10a). 

Given a canonical Hamiltonian operator H{[x,p],t) which satisfies this "Ehrenfest" requirement, the 
substitution of Eq. (10a) into Eq. (9d) permits a series of steps that are analogous to those of the classical 



Eq. (8e), and that likewise result in the explicit algebraic factorization of the divergence operator Vy- 
out of the Eq. (9d) expression for the divergence Vy • j(y,f; [\'ip{t)) , H{[x,p],t)]) of the probabiUty vector 
current density, thereby yielding the probability vector current density j{y,t; [\tp{t)), H{[x,p],t)]) itself as 
the expectation value of the quantization of the classical single-particle scalar density S^^\x — y) times the 
classical particle velocity Vpi?ci([x, p],t). Carrying out this substitution of Eq. (10a) into Eq. (9d), and then 
proceeding in analogy with Eq. (8e), we obtain, 



Vy • j(y, t; [\m),H{[% p], i)]) = -{m\ - y), ^fci([x, p],t)} \m) = 



(Vxc^^^x-y)) ■Vpgc i([x,p],t) m)) = {m\ (Vy^<^ ix-y)) •VpiJei([x,p],t) m)) = 
{m\ Vy [ j(^^x-y)Vpgci([x,p],t)] ^ \m) = 
Vy • [m)\ >i;x-y)Vpi?ci([x,p],t)^ \m)] , 

(lOd) 

where, bccaiise qiiantization and taking the divergence with respect to the non-quantized vector variable 
y are independent linear processes, we can extract the explicitly factored divergence operator Vy- out of 
the quantization — and also, of course, out of the expectation value with the state \ip{t)). Therefore, for 
the subclass of canonical Hamiltonian operators H{[x.,p],t) which satisfy the "Ehrenfest" requirement of 
Eq. (10a), the probability vector current density is given by, 



j{y,t;mt)),H{[%p],t)]) = m)\ <5(^lx-y)Vpi/ei([x,p],t) |^(t)), W 

which is the expectation value of the quantization of the classical single-particle scalar density S^-'^\x — y) 
times the classical particle velocity Vpi?ci([x, p], i). 

Now in addition to those particular "Ehrenfest-subclass" canonical Hamiltonian operators which have 
the form given by Eq. (10c), canonical Hamiltonian operators which have the form, 



H{[%p],t)= p-M^,t) = i(p-vo(x,t)+vo(x,i)-p), (lOf) 

as well satisfy the "Ehrenfest" requirement of Eq. (10a), which we now show in an abbreviated manner, 

[^(^(x - y), i(p • M% t) + M% t) ■ p)] = zn{V-6(^\x - y)) • vo(x, t) = 
in{-VyS^^\Si - y)) • vo(X, t) = -ihVy • [S^^\x - y) Vo(x, t)] , 

and, 



{^(^^x - y), p . vo(x ,t)} = (V,^(^^x-y )) • vo(x, t) = (-Vy^^^^x - y)) • vo(x, t) = 

-Vy [(5(^tx-y)vo(x,<)] =-Vy[S(^\x-y)M%t)]. 

These results make it apparent that for this case as well, the divergence operator Vy- explicitly algebraically 
factors out, paving the way for the resulting probability vector current density to be equal to the expectation 
value of the quantization of the classical single-particle scalar density 5^^\x — y) times the classical particle 
velocity vo(x,t) = VpiJci([x, p], i). 

Combining the special cases of Eqs. (10c) and (lOf), we see that canonical Hamiltonian operators 
H{[x,p],t) which adhere to the "Ehrenfest" requirement of Eq. (lOa) have the form. 



H([%p],t)= i^([p],i) + p-vo(x,t) + y([x],t) = (lOg) 
K{[p],t) + i(p • vo(x,t) + vo(x,i) ■ p) + V{[Si],t). 

For this "Ehrenfest subclass" of the canonical Hamiltonian operators, the divergence operator Vy- always 
explicitly algebraically factors out of the Eq. (9d) expression Vy -j(y, t; [|V'(i))> -'^([x, p]) *)]) for the divergence 
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of the probability vector current density, as we see from Eq. (lOd), and the consequent probabihty vector 

current density j(y, [\ijj{t)), H([Si,p\,t)]) itself is always given by the expectation value of the quantization 
of the classical single-particle scalar density 6^^\x — y) times the classical particle velocity Vpifci([x, p], f), 
as we see from Eq. (lOe). 
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